function path = pathl(vecx,vecy,vecz)
path = sum(sqrt(diff(vecx).^2+diff(vecy).^2+diff(vecz).^2));